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We calculate the 0(a s ) corrections to the production of a hard and isolated photon 
accompanied by one or two jets in deep inelastic lepton nucleon scattering at HERA. 
Numerical results are presented and the potential of this process for studies of parton 
distribution functions is discussed. 
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1 Introduction 

The production of hard photons in hadronic processes is an important testing ground for 
QCD. Since the photon does not take part in the strong interaction, it is a 'direct' probe 
of the hard scattering process. Direct photon production in 7p |]J and in pp collisions 
provides a means to determine the strong coupling constant a s and has been used to extract 
information on the parton distributions, in particular the gluon density in the proton ||. 
In e + e~ annihilation ||], measurements of photon radiation in hadronic Z decays at LEP1 
have provided important independent information on the electroweak couplings of up and 
down quarks to the Z boson [EL ||. Moreover, final states containing a photon are an 
important background for many searches for new physics and a good knowledge of the 
standard model predictions for direct photon production is therefore required. 

At HERA, radiative deep inelastic scattering, ep — > e'yX, with photons collinear to the 
incoming electron has been used to obtain a measurement of the structure function F 2 at 
low values of the momentum transfer Q 2 ]7| . Also the first observation of hard non-collinear 
photons at Q 2 = 0, i.e. in photoproduction has been reported recently ||. With increasing 
luminosity this measurement is expected to contribute information on the parton content 
of the photon and the proton. By contrast, direct photon production at large Q 2 would be 
sensitive to the parton distributions in the proton only. The information obtained this way 
would be complementary to the F 2 measurement from inclusive deep inelastic scattering, 
since up and down quarks contribute with different weights. Typical cross sections for the 
production of hard photons in deep inelastic scattering with Q 2 > 10 GeV 2 are of the order 
of 10 pb. With a luminosity of 50 pb" 1 one thus expects statistical uncertainties of the 
order of 5 % and a measurement of differential cross sections seems feasible. 

Whereas next-to-leading order calculations for direct photon production are available 
for photoproduction |j], H, pp collisions [|10| , as well as for e + e~ annihilation |Tl|, J12] , a 



corresponding calculation for deep inelastic ep scattering was still missing. In this work 
we study the 0(a s ) corrections to the process ep — > e^X at large Q 2 . Since hard photon 
production is a process of relative order a e = 1/137 with respect to the total deep inelastic 
scattering cross section, we expect sizable event rates only at moderately large Q 2 and 
restrict ourselves therefore to pure photon exchange, i.e. Z-exchange contributions are 
neglected. The calculations will be organized in such a way that the hadronic final state 
can be separated into 7 + (1 + l)-jet and 7 + (2 + l)-jet topologies (the remnant being 
counted as "+1" jet, as usual). Our approach is thus analogous to that in calculations of 
(2 + 1)- and (3 + l)-jet cross sections in deep inelastic scattering where a gluon is replaced 
by a photon [|13|]. 7 + (2 + l)-jet events originate through the emission or absorption of a 
gluon. Therefore the ratio of 7 + (2 + l)-jet and 7 + (1 + 1 ) - j et events is sensitive to the 
value of the strong coupling constant a s and to the gluon distribution. 

In addition to perturbative direct production, photons are also produced through the 
'fragmentation' of a hadronic jet into a single photon carrying a large fraction of the 



jet energy [14|. This long-distance process is described in terms of the quark-to-photon 



and gluon-to-photon fragmentation functions. The necessity for taking into account non- 
perturbative contributions is signaled by the presence of singularities showing up in a 
perturbative calculation. These singularities are related to collinear photon-quark configu- 
rations. The factorization theorem of QCD guarantees that all singularities can be absorbed 
into well-defined universal parton-to-photon fragmentation functions, the remainder being 
calculable in perturbation theory. 

In practice, a measurement of direct photon production is feasible only when isola- 
tion conditions are imposed on the observed photon in order to reduce various hadronic 
backgrounds, in particular from two-photon decays of 7r°. The contribution from non- 
perturbative parton-to-photon fragmentation, being related to collinear photon emission 
from partons, can be reduced by isolation requirements, but is not completely removed. 
Again, in a perturbative calculation, this is related to the presence of singularities. In fact, 
if one tries to model the experimental isolation conditions by imposing cutoffs on parton- 
level jets, one can not exclude contributions due to soft quarks having emitted a hard 
collinear photon; the soft quark may appear only as part of a parton-level jet, but not as a 
separate, observable jet which can enter the isolation conditions^]. The implementation of 
photon isolation is particularly non-trivial in a calculation including 0(a s ) contributions 



since the isolation conditions affect the available phase space for gluon emission [19]. As 
a consequence, the parton-to-photon fragmentation functions may have to be modified for 
isolated photon production and higher-order corrections may turn out to be large and to 
require their resummation. 

In the present work we adopt a simpler approach where the fragmentation contributions 
are ignored completely. The photon-quark collinear singularities then have to be removed 
by explicit parton-level cutoffs. The dependence of the final results on these cutoffs (dis- 
cussed in section 4 below) will indicate to what extent the quark-to-photon fragmentation 
function would contribute in a more systematic treatment. 



2 The Leading-Order Process 

In leading order (LO), the production of photons in deep inelastic electron (positron) 
proton scattering is described by the quark (antiquark) subprocess 

e(pi) + q(ps) -* e(p 2 ) + g(p 4 ) + l(p 5 ) (1) 

where we have given the definition of the particle momenta in parentheses. The momentum 
of the incoming quark is a fraction of the proton momentum pp\ p^ = £pp. The proton 



lr The problem is most easily visible in e + e~ — > 7 + 1-jet, where already at leading order photon-jet 
isolation does not remove the photon-quark collinear singularity |12j| . A next-to- leading order calculation 
|L5j shows features typical for a next-to-next-to-leading order calculation. Measurements of the quark-to- 
photon fragmentation function in e + e~ — > 7 + 1-jet had been proposed in Refs. |l6|, and were described 
in Refs. 111. 



remnant r carries the momentum 

Pr = (1 - OPP (2) 

and hadronizes into the remnant jet so that the process (0) gives rise to 7 + (1 + l)-jet 
final states. The momentum of the hadronic final state, i.e. the (1 + l)-jet system, is 
Pp + Pi — Pi — Pn an d its invariant mass W is given by 

W 2 = (pp+ Pl -p 2 -p 5 ) 2 . (3) 

We will use the well-known kinematic variables for deep inelastic scattering 

Q 2 = ~(Pi -P2) 2 , x=- — -, r, y= — , s = (p 1 +p P ) 2 } (4) 

2pp{Pi-P2) xs 

determined by the momentum of the scattered lepton. Because of the presence of the 
photon in the final state, large Q 2 does not guarantee large W and we will have to require 
explicitly W > W m in in order to stay in the deep inelastic regime where a perturbative 
treatment can be expected to work. Apart from this, we will also apply cuts on the variables 
x, y and Q 2 since we ask for an observable scattered electron. These latter cuts remove 
direct photon production in photoproduction. 

Both leptons and quarks emit photons. The subset of Feynman diagrams where the 
photon is emitted from the lepton ("leptonic radiation") is gauge invariant and can be 
treated separately. Similarly, the Feynman diagrams with a photon emitted from the quark 
line is called "quarkonic radiation". There is also a contribution from the interference of 
these two parts. For tests of QCD the interest is in those contributions where the photon 
is emitted from quarks and leptonic radiation is viewed as a background. 

Radiative deep inelastic scattering appears as a contribution to QED radiative correc- 



tions (see for example Ref. [20] and references therein). In this case the emitted photon 
remains undetected and singularities due to soft and collinear photons have to be canceled 
by taking into account virtual 0(a) corrections to non-radiative scattering eq —>■ eq. Here 
we are interested in events with an observable photon, i.e. we restrict ourselves to the case 
where the energy of the photon E^ = E 5 is sufficiently large, 

Also, the photon should be spatially separated from all other particles: 

#7,i > #sep, (6) 

where # 7i j is the angle between the momenta of the photon and particle i (= 1,2,3,4 for 
the leading-order process (HD and similarly for the next-to-leading order processes specified 
in section 3 below). In particular, the photon is not allowed to be emitted close to the 
beams: 

"min ^ "7 *^. "max • \' ) 



These cuts remove all photonic infrared and collinear singularities. Instead of using the 
angle 9 y> i, photon separation from final state particles can also be imposed by cuts on the 
invariant masses 
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or, normalized to the invariant mass of the hadronic final state, 

■* = £■ 0» 

The condition 

y^ > yo (io) 

(i = 2, 4, r) is more comfortable for the analytic calculation, but less suited to experimental 
requirements. Since we will perform the phase space integration with the help of Monte 
Carlo techniques, we are not restricted to one specific choice of isolation criteria, but we 
can apply a combination of the above cuts as will be described below. 

At lowest order, each parton is identified with a jet and photon-parton isolation corre- 
sponds to the isolation of the photon from an observable jet. With isolation cuts, parton- 
to-photon fragmentation does not contribute at this order. 



3 0(ol s ) Corrections 

At next-to-leading order (NLO), processes with an additional gluon, either emitted into 
the final state or as incoming parton, have to be taken into account: 

e(pi) + q(p 3 ) -> e(p 2 ) + q(p 4 ) + j(p 5 ) + gfa), (11) 

e(Pi) + g(pa) -* e(p 2 ) + g(p 4 ) + t(pb) + QiPe), (12) 

where the definition of momenta is again shown in parentheses (see Fig. [J]). In addition, 
virtual corrections (one- loop diagrams at 0(a s )) to the process ([!]) have to be included. 

The amplitude for purely leptonic radiation at order 0(a s ) factorizes into a leptonic 
tensor for e — > cyj* and a hadronic tensor including next-to-leading order QCD cor- 
rections. Both parts are well-known and their combined contribution to deep inelastic 
scattering is included for example in the Monte Carlo program DJANG06 ||21|| . For the 



0{a s ) corrections to quarkonic radiation and in particular the lepton-quark interference, 
a representation in terms of a leptonic and a hadronic tensor is not suitable. The corre- 
sponding complete matrix elements including the leptonic and hadronic vertex have been 
obtained with the help of form [K3] and are given in 



Whereas the LO process leads to the appearance of events with a photon and one 
current jet, 7 + (1 + l)-jets, in higher orders additional jets can be produced: the processes 
( |TT| , |T2"D contribute both to the 7 + (1 + l)-jetcross section, as well as to the cross section for 





Figure 1: Examples of Feynman diagrams for eq — > eqyy (a) and eg — > eqq^ (b) with the 
definition of momenta. 

7 + (2 + l)-jets, depending on whether the quark-gluon or quark-antiquark pair in the final 
state appears as one single jet or as two separated jets. The two cases can be identified by 
comparing the scaled invariant masses of parton pairs with a jet resolution parameter y J : 
two partons (i,j) with i,j = 4, 6, r are supposed to lead to 2 jets if 

vn > y J - (is) 

Also the remnant r is treated as a parton and a quark, antiquark, or gluon in the final state 
is recombined with the remnant into one jet if yi r = -jrya is smaller than y$ . Similarly, 
photon isolation can be imposed with the help of cuts on the scaled invariant masses. 

In the phase space region where two jets can not be separated, the matrix elements 
become singular. These singularities appear when one of the partons becomes soft or when 
two partons become collinear to each other. The singularities can be assigned either to 
the initial state or to the final state (ISR: initial-state radiation, FSR: final-state radia- 
tion). The FSR singularities cancel against singularities from virtual corrections to the 
lower-order process. For the ISR singularities, this cancellation is incomplete and the re- 
maining singular contributions have to be factorized and absorbed into renormalized parton 
distribution functions 



To accomplish this procedure, the singularities have to be isolated in an analytic cal- 
culation, e.g. with the help of dimensional regularization. The application of dimensional 
regularization is, however, not feasible for the complete cross section of the higher-order 
processes. Therefore we use the so-called phase-space slicing method |25) to separate those 
regions in the 4-particle phase space which give rise to singular contributions. A separation 
cut yl is applied to the scaled invariant masses y^ and chosen small enough, such that 
the calculation can be simplified by neglecting terms of the order 0[y'l\ Contributions 

6 



from phase space regions where one of the y^ is smaller than y$ are singular and have to 
be combined with the one- loop corrections to obtain a finite result. The sum of these two 
contributions defines the cross section for events where two partons are recombined into a 
parton-level jet (parton-level (1 + l)-jet events). The contributions where all y^ are bigger 
than i/q are related to final states with three separate partons (parton-level (2 + l)-jet 
events). The latter are free of singularities and can be calculated with the help of Monte 
Carlo techniques. 

As known from similar calculations (e.g., for the (non-radiative) jet cross sections in 
DIS p6fl), the phase space slicing parameter y$ has to be chosen very small, of the order 



of 10~ 3 or smaller, in order to allow for the neglect of terms of order 0(?/q). Therefore, y$ 
can not be identified with the y-cut of a jet algorithm applied in an experimental analysis. 
There, due to experimental restrictions, y cannot be reduced to values below O(10~ 3 ). 
In addition, a fixed-order calculation may give unphysical, i.e. negative (1 + l)-jet cross 
sections for too small values of y (see the curves labeled with S in Figs. @a and ^|a below). 
The Monte Carlo approach, however, allows to apply a jet algorithm to the parton-level 
events, i.e. to recombine 2 partons in the parton-level (2 + l)-jet events according to a jet 
algorithm using y-cuts y J for the separation of jet pairs (similarly: y 1 for the separation 
of a jet and a photon) with values as appropriate for the given experimental situation. 

The calculation thus proceeds through two subsequent steps: First, phase space slicing 
is applied with a small y-cut y$ of the order of ^ 10~ 3 to accomplish the cancellation of 
singularities. This step relies on analytic calculations. Secondly, a jet algorithm is applied 
with experimentally realizable, i.e. large enough values y J and y 1 of the order of 0.01 — 0.1. 
The second step is performed during the Monte Carlo integration. 

The singular contributions for the process eq — > eqy-f involve the following factors 

(14) 




those for eg — > eqq^y contain the factors 



• 



2/36 2/34 



(15) 



Terms containing l/t/3j, i.e. the momentum j?3 of the incoming parton, are associated 
to initial-state singularities, terms that do not, to final-state singularities. Contributions 
involving the product of an ISR and an FSR factor, as for example the factor l/j/36j/46, can 
be separated by partial fractioning, 

11111 (16) 



UijVik yij Vij + Vik Vik ytj + ytk 



so that all singular contributions can be associated either to the initial state or to the final 
state. Note that the denominator y^ + y^ introduced by partial fractioning can become 



zero only if both y^ = and i/ik = at the same time; since configurations where all three 
partons 3, 4, and 6 are collinear with each other are excluded by the cut on W, this is 
possible only for p t = 0. Therefore, for a contribution containing the pole factor 1/yij, we 
can separate the phase space into three regions: 

• yij < y'l- This region contains the infrared singularity at y^ = yu~ = 0, as well 
as the collinear singularity at y. t j = 0, y^ > and leads to singular contributions, 
i.e. 1/e and 1/e 2 poles in dimensional regularization. The double-poles 1/e 2 and 
parts of the single-poles 1/e cancel with corresponding singular contributions from 
virtual corrections. The remaining 1/e-pole contributions are associated to the initial 
state, can be factorized, and are absorbed by renormalizing the parton distribution 
functions. The analytical integration over this phase space region is performed with 
the approximation of small j/q, i.e. neglecting terms of 0(?/q). This contribution will 
be denoted by "S" (singular) below. 

• Vij > Vo an d Vik > Vq with only parton-level 7 + (2 + l)-jet events, denoted by "R" 
(real corrections); 

• Vij > Vo and y^ < y$. Here, the result is non-singular (therefore denoted by "F", 
finite) but does not vanish with y£ — ► 0, contrary to naive expectations. Its contri- 
bution is calculated numerically. It is non-negligible in particular for terms related 
to ISR singularities. 

The integrals needed for the singular contributions are written in a Lorentz-invariant 
form as tensor integrals which can be reduced to a few basic scalar integrals with the help of 
analytic programs like mathematica or form. More details are given in |23|. The remaining 
phase space integrations are performed with the help of Monte Carlo techniques. The three 
contributions S, R, and F are treated separately, each with appropriate mappings of the 
respective integration variables to improve the numerical stability of the calculation. 

As discussed in the introduction, in the present work we do not factorize and subtract 
those photon-parton collinear singularities which have to be absorbed into parton-to-photon 
fragmentation functions. Instead, we remove all singular contributions by keeping isolation 
cuts at the parton level. As stated in the introduction, care has to be taken that the 
isolation criteria do not restrict the phase space for gluon emission since this would destroy 
the cancellation of singular contributions. Therefore, in the first step of the calculation 
described above, we require the photon to be isolated from the quark (antiquark) by the 
cut 

Vu > Vo ( 17 ) 

with % = 3,4 for eq — > eqyy and % = 4, 6 for eg — > eqcpf. The cut is not applied to 
photon-gluon pairs which is possible since gluons do not emit photons and there is no 
singularity related to y gi . This definition of photon isolation at the parton level introduces 
an unphysical parameter (?/q). The sensitivity to yj can be reduced by applying, in the 



second step of the calculation, photon isolation with respect to jets described by cutoff 
parameters which can be used in the same way in the experimental analysis. In order to 
have some freedom when modeling these physical isolation criteria we choose a small value 
for %)q. The dependence on y^ will be discussed below. 



4 Numerical Results 

The results discussed in the following are obtained for energies and cuts appropriate for 
the HERA experiments: the energies of the incoming electron (positron) and proton are 
E e = 27.5 GeV, E P = 820 GeV and 

Q 2 > 10 GeV 2 , W > 10 GeV, 
0.001 < x < 0.5, 0.05 <y< 0.99, (18) 

pT > 5 GeV, 90° < Q~ t < 170°, 6 ie > 10°. 

Note that the emission angle of the photon, # 7 , measured with respect to the incoming 
electron in the HERA laboratory frame, is restricted to the hemisphere # 7 > 90° since 
photon production with # 7 < 90° is dominated by 'uninteresting' leptonic radiation. The 
parton distribution functions are taken from Ref. p8| (MRS (A)). 



The events generated during Monte Carlo integration are 75, 'yqg or ^qq events. A 
simple event analysis is applied to obtain 7 + (1 + l)-jet and 7 + (2 + l)-jet event samples. 
The event analysis consists of two parts: the first part serves to identify the number of 
jets according to a conventional jet algorithm; the second part treats photon isolation. For 
simplicity we choose a jet definition using the normalized invariant masses yij. Since for 
small £ > x the momentum of the remnant and thus yi r can be large even for partons with 
small transverse momentum, we first remove low-p T partons before recombining partons 
to jets. Explicitly we apply the following conditions: 

(1) A final state parton (quark or gluon) is recombined with the remnant if its transverse 
momentum is below a cutoff: 

pj < P l in = lGeV, z = 4,6. (19) 

(2) Two partons are recombined into one jet if 

Vij<y J for i,j = A,6,r (20) 

and all quarks, antiquarks and gluons as well as the proton remnant are taken into 
account when forming jets. If several pairs of partons have y^ below y J , the pair 
with the smallest y^ is recombined first. Several prescriptions for the recombination 



are possible: the energy and 3-momentum of a jet (ij) obtained from pairing partons 
i and j can be obtained by 

E i:i = a{E { + Ej), pij = f3{pi + pj). (21) 

For example in the E-scheme one chooses a = (3 = 1; in the P-scheme one has 
a — Wi + Pj\ /(Ei + Ej), (3 = 1 instead. Since the present calculation is of first order 
in a s , the recombination has not to be iterated. However, the different recombination 
prescriptions become relevant when photon isolation with respect to jets is imposed. 
In addition, the cut on low-pr partons or parton-pairs Eq. (|H]) is affected if a recom- 
bination prescription with j3 ^ 1 is used. Our numerical results will be given for the 
P-scheme. 

(3) Finally, an event is accepted only if the photon is separated from the jets or if the 
photon is accompanied by hadronic energy less than a specified amount, i.e. we 
exclude events with 

y^j < y 1 and Ej > e (Ej + E 7 ) (22) 

where j denotes any jet (i.e., parton or pair of partons) remaining after steps (1) and 
(2) of the event analysis. 

We keep the possibility to use different values for the y-cuts applied to purely hadronic jets 
and to jets containing the photon. In practice, y J and y 1 are taken equal with a typical 
value 0.03. For the photon isolation parameter e we will take the value 0.1 as used in 
experimental analyses ||. Apart from being experimentally unrealistic, the value e = is 
theoretically not allowed since Eq. (|22| ) with e = would restrict the phase space for soft 
partons and consequently destroy the cancellation of corresponding singularities. 

We start with demonstrating the consistency of our approach by showing the depen- 
dence of the 7 + (1 + l)-jet cross section on the phase space slicing cut y$ . Figures ^| and 
|3| show the dependence of the total and partial cross sections for 7 + (1 + l)-jet events. 
For g(g)-initiated processes, the separate contributions S and R depend on log 2 y$ (see 
Fig. ^|a) and the finite contribution F is not negligible. In this case the sum is numerically 
stable in the range 10~ 5 ^ y$ ^ 10~ 3 ; for smaller values the numerical precision decreases 
and for larger values the error from neglected terms of order 0(y^) is not negligible. The 
calculation of (^-initiated contributions can be performed with much smaller uncertainties 
and for much smaller values of y$ , as seen in Fig. |3] since here the dependence on logy^ is 
only linear. Also, the finite contribution is negligible for diagrams with incoming gluons. 
The dependence on y^ at large values above ~ 10 -3 is slightly stronger in this case than 
for g(g)-initiated contributions since terms of order 0(y^) are relatively more important. 
In the following we fix y$ at the value 10~ 4 . 

Figure § shows the residual dependence of the cross sections for 7 + (1 + l)-jets on 
the parton- level photon isolation cut y$, separately for (anti) quark and gluon- initiated 
processes. The y$ -dependence is weak for the case with incoming (anti) quarks showing that 
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Figure 2: Dependence of the 7 + (1 + l)-jet cross section on the phase space slicing cut y^ 
for incoming quarks and antiquarks (y$ = 10^ 4 , y 1 = y J = 0.03J. (a) shows the separate 
contributions and the sum = R + S + F + LO. (b) shows the sum of all contributions, 
including the leading-order cross section, on a larger scale. 
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Figure 3: Dependence of the 7 + (1 + l)-jet cross section on the phase space slicing cut y$ 
for incoming gluons (yl = 10~ 4 7 y 1 = y J = 0.03,). (a) shows the separate contributions and 
the sum = R + S. (b) shows the sum of all contributions on a larger scale. 
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Figure 4: Dependence on the infrared cutoff parameter y$ of the j+(l + l)-jet cross section 
for incoming (anti)quarks (a) and gluons (b) (yl = 10 -4 , y 1 = y J = 0.03/ 



the isolation criteria efficiently reduce the sensitivity to the phase space region where the 
non-perturbative parton-to-photon fragmentation functions would contribute. For the g- 
initiated processes, the sensitivity to y$ is larger. In this case all final-state partons (q and q) 
can emit a photon. Since the isolation condition is applied to jets, the singularity associated 
to configurations with soft (anti) quarks having emitted a hard collinear photon is removed 
only with the help of the parton-level cut Eq. (|TTD. For processes with incoming quarks, 
only a small subset of diagrams leads to singularities for similar configurations. We choose 
y$ = 10 -4 in the following. This value is small enough compared with experimentally 
realistic values for y" f ^ O(10~ 2 ) so that contributions where a quark or an antiquark 
determines the momentum of a jet, become insensitive to y$. Also, much larger values 
would lead to an unphysical negative cross section for the ^-initiated subprocess. In a 
more systematic treatment, the ^"dependent terms in our calculation would be replaced 
by contributions from parton-to-photon fragmentation functions. In our present approach, 
however, the unwanted dependence on y$ has to be viewed as an unavoidable source of a 
theoretical uncertainty. The 7 + (1 + l)-jetcross section at Q 2 ^ 100 GeV 2 is affected by 
this at the level of 20 % (see Fig. || below). At larger Q 2 , the influence of the y$ -dependent 
gluon-initiated contribution is reducedQ 

In Fig. ^ we show the differential cross section da/dd^ (sum of 7 + (1 + l)-jets and 
7 + (2 + l)-jets) in the range 10° < 6* 7 < 175°. Apart from the extended range of photon 

2 This can be compared with the case of e + e~ — v 7+hadrons where the total cross section has little 
sensitivity to the parton-level photon isolation cut for not too large y, but the 7 + 1-jet rate has a non- 
negligible dependence on j/q fl27|| . 
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Figure 5: Differential cross section da/d9 1 forep — > e7 + (1 + 1) -jets and ep — > e7 + (2 + l) 
jets with y J — y 7 — 0.03. Full histogram: lowest order, dashed histogram: including 0(a s 
corrections. Cuts are explained in the text. 



emission angles all cuts given in ( lq) are applied. The majority of photons is produced 
with small angles, i. e. close to the direction of the incoming lepton. For leptonic radiation, 
QCD corrections reduce the cross section by ~ 10 % for the phase space region under 
consideration. By contrast, at large emission angles, dominated by quarkonic radiation, 
the cross section receives positive QCD corrections. In the following we restrict ourselves 
again to this "signal" region 9 1 > 90°, i.e. the proton hemisphere in the HERA laboratory 
system. 

The Q 2 -dependence in this restricted phase space region is shown in Fig. |6|. The cross 
section is shown separately for g(g)-initiated and gluon-initiated contributions giving rise to 
7 + (l + l)-jet and 7+(2+l)-jet events using y J = y 1 = 0.03. The 7 + (l + l)-jet contribution 
is dominant for these y-cut values with a maximum in the lower Q 2 range, whereas the 
cross section for 7 + (2 + l)-jet events is flatter and extends to larger Q 2 . Incoming 
gluons contribute only roughly 10 % to the total cross section. Since the distributions for 
incoming quarks and incoming gluons are not very different, it seems difficult to utilize 
radiative deep inelastic scattering for a measurement of the gluon distribution. We also 
checked that using other parametrizations of parton distribution functions (like those of 
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Figure 6: Differential cross section Q 2 do~/dQ 2 for ep — > cy + (n + l)-jets with y J = y 1 = 
0.03. Upper full histogram: lowest order, dashed histogram: contribution from incoming 
(anti)quarks to 7 + (1 + l)-jets, dotted histogram: contribution from incoming (anti) quarks 
to 7 + (2 + l)-jets, dash-dotted histogram: incoming gluons for 7 + (1 + l)-jets, lower full 
histogram: incoming gluons for 7 + (2 + l)-jets. Photons are restricted to 90° < Q 1 < 170° 
and other cuts are explained in the text. 



Refs. [£9]) do not lead to significantly different shapes of distributions. Only the total cross 
sections vary by 0(1 — 15 %). 



The rate of 7 + (2 + l)-jet events, 



R 



7,2+1 



a( 7 + (2 + 1) -jets) 



(7(7 + (1 + l)-jets) + o-(7 + (2 + l)-jets) 



(23) 



increases towards smaller values of y J (see Fig. [F|) and becomes equal to the 7 + (1 + l)-jet 
3 , the precise value depending on y~> '. The dependence on y 1 is weaker; in 



rate at y J ^ 10 



particular, for y 1 ^ 0.02 the 7 + (2 + l)-jet rate is almost independent on y 1 . The reduction 
of the cross sections with increasing y 1 is stronger for 7 + (2 + l)-jet events at large values 
of y J than at small y J , relative to the 7 + (1 + l)-jet cross section, i.e. the ratio -R 7 ,2+i 
increases with increasing y 1 at large y J whereas it decreases with increasing y 1 at small y J . 
For completeness we present the y-cut dependence of the 7 + (1 + l)-jet and 7 + (2 + l)-jet 
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Figure 7: 7 + (2 + I) -jet rate -R 7 ,2+i (is a function of the jet cut y J for three different values 
of the photon isolation cut y 1 . 
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Figure 8: 7 + (1 + l)-jet and 7 + (2 + l)-jet cross sections as a function of jet cut and 
photon isolation. 
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cross sections in Fig. |[ Note that also the total cross section has a dependence on y J , as 
can be seen from the sum of the results shown in Figs. ||a and b. The jet algorithm not 
only defines the classification of the hadronic final state into 7 + (1 + l)-jet or 7 + (2 + l)-jet 
events, but also affects the overall phase space boundaries: smaller values of y J allow the 
jets to be closer to the remnant jet so that the cut against \ow-p T partons has a stronger 
effect. 




Figure 9: Dependence on the renormalization and factorization scales /if, = //|, = 
the total cross section (i.e., the sum of ^+(1+1) -jets and ^+(2+1) -jets) fory J = y^ 



fQ 2 of 

= 0.03. 



The leading order cross section depends on a factorization scale fip via the scale entering 
the parton distribution functions qi(x,fi 2 F ). At next-to-leading order, there is an explicit 
scale dependence in the 7 + (l + l)-jet cross section through factorization of the initial state 
singularities which partly compensates the scale dependence from the parton distribution 
functions. In addition, the explicit factor a s depends on the renormalization scale. For 
simplicity we identify the two scales, which could in principle be chosen independently 
from each other. Figure |9] shows the scale dependence of the leading-order and the next- 
to-leading order cross sections where we have used //§, = fjL 2 R = fQ 2 . Varying / between 0.25 
and 4, the total cross section for ep — > e^X shows good stability within a few percent. Also 
from this figure one can infer that the 0(a s ) corrections vary between 20 and 30 %. A 5 % 
measurement of the cross section would therefore correspond to a 20 to 30 % measurement 
of a s , assuming negligible uncertainties from parton distribution functions. 
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Radiative ep scattering is complementary to usual deep inelastic scattering since up 
and down-type quarks contribute with different weights to the cross sections. In the usual 
structure function F 2 , the sums of u- and d-type quarks, U = u + c + u + c and D = 
d+s+b+d+s+b enter with the relative factors e\ : e\ = 4 : 1 whereas for the 
contribution from quarkonic radiation to ep —>■ ejX this ratio is e* : e\ — 16 : 1. In 
principle, a common analysis of non-radiative and radiative scattering performed with 
high enough precision, would allow a determination of U and D separately. We therefore 
investigated the dependence of the total cross section for ep — > e^X on the ratio U/D. 
In order to keep the well-constrained structure function F 2 unchanged, we modified the 
parton distributions by the following prescription: 




(24) 
AiU/D))' 

By this the combination e A JJ + e\D is replaced with [1 + 3(1 - 8 d )/(l + 1QU/D)} x (e A u U + 
e d D). With U/D ~ 1.5, a typical value at x ~ 0.1, one expects a 12% reduction for 
Sd = 2 and a 6 % enhancement for 5d = 0.5. In fact, the true change of the cross section is 
smaller (—5.7 % and +2.3 % with the cuts (H)) since additional contributions from leptonic 
radiation and quark-lepton interference, which are not proportional to the fourth power of 
the quark charges, are not negligible even in the 'signal' region # 7 > 90°. It thus seems 
unlikely that with respect to a determination of U/D radiative deep inelastic scattering 
could become competitive with classical analyses like that of the difference of proton and 
neutron cross sections or the W charge asymmetry in pp — > W ± + X. 



5 Summary 

We have described a first next-to-leading order calculation of isolated photon production 
in ep scattering at large Q 2 . Apart from providing a sound basis for testing QCD in direct 
photon production, our results improve the knowledge of standard model predictions as a 
source of background for searches for new physics. We have discussed numerical results for 
7 + (1 + l)-jet and 7 + (2 + l)-jet cross sections at HERA. Corresponding measurements 
will provide valuable information that will allow to further constrain parton distribution 
functions, in particular when combined with results from other experiments. It still has 
to be investigated which kinematical variable would be best suited to obtain the highest 
sensitivity on the gluon distribution, the U/D ratio, or the strong coupling constant a s . 
For example, in photoproduction the distributions with respect to the photon rapidity or 
photon transverse momentum (in the HERA laboratory or in the 7*p center-of-mass frame) 
turned out to be good choices. One should also expect that the theoretical uncertainties 
due to the parton-level cutoff yl could be further reduced by optimizing the analysis with 
respect to kinematical cuts, the jet algorithm {e.g., a cone algorithm as used in the study 
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of (2 + 1)- or (3 + l)-jet events at HERA p0[ ) and modified photon isolation prescriptions 



(e.g., the so-called "democratic" clustering procedure fig) ). 
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